Sign-time distributions for interface growth 
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We apply the recently introduced distribution of sign-times (DST) to non-equilibrium interface 
growth dynamics. We are able to treat within a unified picture the persistence properties of a 
large class of relaxational and noisy linear growth processes, and prove the existence of a non-trivial 
scaling relation. A new critical dimension is found, relating to the persistence properties of these 
systems. We also illustrate, by means of numerical simulations, the different types of DST to be 
expected in both linear and non-linear growth mechanisms. 
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The notion of persistence, or the statistics of first 
passage events, has been a powerful conceptual tool in 
studying stochastic non-Markovian processes in many re- 
search areas of physics, engineering, statistics, and ap- 
plied mathematics. In this Letter we apply a persistence- 
related concept, the distribution of sign-times, or DST, 
(defined below) to the problem of kinetic surface rough- 
ening in non-equilibrium interface growth dynamics 0. 
We believe that the ideas, described in this paper, could 
become an extremely useful conceptual and practical 
tool in characterizing surface growth dynamics, rival- 
ing the dynamic scaling ideas currently used in study- 
ing kinetic surface roughening. Depending on the spe- 
cific issues of interest, our proposed DST technique may 
actually be more powerful and informative than the 
currently fashionable dynamical/roughness/growth ex- 
ponent based characterization of dynamical surface mor- 
phologies. 

One of the main themes in the theory of non- 
equilibrium interfaces is grouping the interface roughen- 
ing phenomena within 'universality classes'. This classi- 
fication of scenarios is based on calculating the dynamic 
scaling properties of the surface correlation function [jj] . 
On the other hand, in non-equilibrium interface growth 
experiments, one might also be interested in morphology 
stability issues which can in fact be formulated as first 
passage type questions: what is the probability that a 
mound (or crevice) will survive as a mound (crevice) for 
a given period of time tl How does this probability decay 
in time, etc.? These type of questions however are not 
simply delineated by such a correlation function. 

Another open theoretical problem is to establish a cor- 
respondence between discrete solid-on-solid (SOS) mod- 
els and continuum Langevin equations beyond the equal- 
ity of exponents. For example, based on structure factor 
measurements, the authors in Ref. j| claim that the SOS 
model they introduced does not only belong to the same 
universality class as the noisy Mullins equation but it 
is described exactly by it. Our approach proposed in the 
present Letter (which is not based on direct measurement 
of the correlation function) supports that claim. 



It would be useful therefore, to study statistical quan- 
tities that are directly sensitive to the structural and 
morphological properties of interfaces (e.g., formation of 
mounds) and to the dynamics of these structures (e.g., 
coarsening). In this Letter, we propose that such infor- 
mation may be inferred from the DST, which has recently 
been introduced in the context of the persistence proper- 
ties of simple coarsening systems and the diffusion equa- 
tion First passage time or persistence problems 
have been the focus of intensive research for the past 
few years, producing a series of analytic and numerical 
results with applications to the Ising and Potts models 
||, the diffusion equation Q, phase ordering 0, inter- 
face kinetics [§), etc. and experiments on liquid crystals 
and soap froths (see references in The central issue 
of persistence concerns the probability of an event never 
occurring (up to a certain time t). It is very restrictive 
by definition, and good statistics from numerics or exper- 
iments may be extremely hard to obtain. The recently 
introduced ||[|| DST is practically more accessible, and 
as a limiting case produces the persistence probability. 

The DST is essentially a histogram performed on the 
sign of the fluctuations and simply measures the prob- 
ability of the fluctuations having been in the positive 
domain for a total time r in the given time t of the pro- 
cess. Obviously for r = t we obtain the usual persistence 
probability, which we denote by P+(f), and for r = we 
obtain the probability of the fluctuations having never 
been in the positive domain, i.e, to have been always in 
the negative domain, P_ (t) . The distinction between the 
persistence of fluctuations in the positive domain and in 
the negative domain becomes important in the case of 
nonlinear models ||. We shall refer to these as 'positive' 
and 'negative' persistence, respectively. 

The sign-time for an interface on a (i-dimensional sub- 
strate is the stochastic variable defined by: 



T(x,t) = / dt'H(h(x,t)) , 
o 



(1) 



where H is the Heaviside step function and h(x, t) is the 
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height of the interface measured with respect to the av- 
erage height. Since h is a random variable (due to its 
coupling to the noise) the sign-time will be described by 
a probability distribution - the DST. For a system with 
translation invariance, the statistics of t will not depend 
on the location x, and so the DST may be written as 

S(r,t) = (5(r-T(0,t))) , (2) 

where (.) indicates the average over the noise. Some 
properties of S are: i) it is defined on the interval 
< r/t < 1; ii) for interface growth with h —* —h sym- 
metry, S will be symmetric about rjt = 1/2; iii) the 
tails of the distribution give the persistence probabili- 
ties: P+(t) = /* dr S(t, t) and P_(i) = J Q e dr S(t - r,t), 
where e <C t is a microscopic time scale (of the order 
of the fastest temporal scale in the interface dynamics 
). These probabilities are expected to have a power law 
decay, defining the corresponding persistence exponents 
0±: P±(t) ~ [l/t){e/t) 0± - x , iv) the shape of S contains 
information about whether the growth is rare event dom- 
inated or not. 

In the spirit of ref. || we first consider the following 
class of stochastic linear equations: 

d t h = -v(-v 2 y/ 2 h + t , (3) 

with flat (h(x, 0) = 0) initial condition, where £ is a gen- 
eral noise term which may represent the 'pure determinis- 
tic' case via the choice £(x, t) — 5(t)r)(x, t) or the regular 
'noisy' case with £(x, t) = rj(x,t) where rj is a Gaussian- 
distributed noise possibly with spatial correlations. Wc 
consider the following three choices for rj: 1) white noise 
with correlator (r)(x,t)r}(x' ,t')) = 2DS(x-x')S(t-t'), 2) 
volume conserving noise (rj(x,t)r](x' ) t')) = — 2DV 2 5(x — 
x')S(t — t') and 3) long range spatially correlated noise 
(r)(x,t)r)(x',t')) = 2D\x-x'\P- d S(t-t'), p < d. For 
example, the Edwards- Wilkinson (EW) model may be 
recovered by setting z = 2 in Eq. and by applying 
white noise; likewise the noisy Mullins equation corre- 
sponds to setting z = 4 J3j. We write Eq.(Q) through the 
higher moments of DST as 

oo 

#W) = £ / ^t^)l {[T ^ Mrh (4) 

n =°-cc 

where we have introduced a frequency representation of 
the delta function, and expanded in powers of the sign- 
time TJ. We shall enter into no technical details here 
on how to proceed with calculating the moments of the 
DST. We present only the final form that we obtained for 
the n th order moment normalized by t n (p n = {(T~/t) n )): 

n 1 00 

7~\ J 27r J e k + ia k 

K — 1 — oo 



where the limits of e k — > + are to be taken, and 

!{x + 2/)~ 7 , deterministic case 
min r x,v) a , ^ o W-r • ( 6 ) 
J du{x + y — 2u) 1 , noisy case 
o 

with < x, y < 1, and 7 being given by 1) 7 = d/z 
for the deterministic case and for white noise, 2) 7 = 
(d+2) I z for volume conserving noise, and 3)7= (d—p) / z 
for long-range correlated noise. We make the following 
observations from Eqs. (Q)-@- First, the DST obeys the 
exact scaling form 

S i d z \r,t) = - t F 1 Q, 0<r<i, (7) 

for all values oft (//„ is t- independent) . Second, the 'ma- 
terial parameters' v and D do not appear in the DST. 
Third, the three numbers (d, z, p) appear in the DST (for 
any t) only through their combination in 7 = j(d, z, p). 
Thus, the persistence exponents (which are contained 
within the DST) will also only depend on d, z and p 
through the exponent 7. [This appears to be implicitly 
understood in ref. Q where persistence is measured as a 
function of the growth exponent (3 = max(0, (1 — 7)/2)]. 
A similar scaling property for the persistence exponents 
is also true for the deterministic case. 
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FIG. 1. DST's {F 1 (x) = tS(xt,t)) for a) {1, 2} (thick line) 
at t = 0.25 x 2 • 10 3 obtained on a grid of L = 2048 sites and 
averaged over 2 • 10 3 runs; {2, 4} at t = 0.01 x 256 (diamonds) , 
and t = 0.01 x 4096 (thin line), on a grid of 1024 x 1024 
shown for a single run, b) the SOS large curvature model 
(dots) on a lattice of L — 10 4 at 8192 steps, averaged over 
100 runs; and for {1, 4} measured on a grid of L — 2048 sites 
at t = 0.05 x 2 • 10 3 , and averaged over 2 ■ 10 3 runs. 

For simplicity of the notation, instead of (r, t) (and 
6^) we will use S y (T,t) (and 7 ). Let us consider as 
an example the generic case of white noise, for which 
7 = d/z. According to the above, for any model for 
which, e.g., d/z = 0.5, the DST (and thus the persis- 
tence properties) will be identical to that for the EW 
model in one dimension. We compared the numerically 
obtained DST's for {1,2} (meaning d = 1, z = 2) and 
{2, 4} . According to the above, one should observe iden- 
tical DST's. The {1, 2} DST was measured using a stan- 
dard discretization scheme, see Fig. la). The numerical 
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integration of the case {2, 4} is less straightforward. We 
used the simplest discrete scheme for modeling the op- 
erator V 4 in d = 2, as more sophisticated schemes were 
actually less stable under the influence of additive white 
noise. A very small integration time step of St = 0.01 was 
used to ensure stability. We observed a long transient in 
which the DST was actually concave, in contrast to the 
convex DST for the case {1, 2} (d = 1 EW model). After 
10 3 iterations the DST began to turn over, and even- 
tually settled into a convex shape, closely matching the 
{1,2} DST. This illustrates the sensitivity of the DST 
to lattice effects, which may be a very useful property 
if one is investigating physics which is itself sensitive to 
the underlying lattice. In fig. lb) we show the DST's ob- 
tained from numerical integration of the case {1,4} (the 
d = 1 noisy Mullins equation) and numerical simulation 
of the SOS large curvature model These two models 
are expected to have very similar properties, and indeed 
their DST's are almost indistinguishable. 

It follows from Eqs. (jj)-© that the knowledge of the 
second moment /12 uniquely determines the value of 7, 
and therefore also z (in a given dimensionality). It is 
possible to evaluate analytically the second moment \xi- 
We find \i2 — 1/2 — G (7), where for the deterministic 
case 



SJ T ,t) = -s(l-- 

jy ' t \2 t 



G( 7 ) = f /da^ 
w ; 4vr / 1 + a 



and for the noisy case: 



1 + a 
2^72 



2- 



I -1/2 



- 1 



(8) 



Gh) = -7- I da arctg, / ^ 



- 1 



(9) 



The second moment in both cases is a monotonic func- 
tion of 7 and therefore the knowledge of one determines 
the other; a property useful in deciding whether a mea- 
sured DST can indeed be described by a process like (|^). 
For example, one may obtain from numerical or experi- 
mental measurements a symmetric DST, from which one 
may compute fj,2- One can test therefore if the process 
generating the measured DST can be described by Eq. 
(||): one determines z using the above procedure, and 
then simulates Eq. (|^) with the corresponding value of 
7, thus generating a new DST. If the two DST's are very 
close or coincide then the assumption that the physical 
process may be modeled by Eq.(||) is valid, just as in the 
SOS large curvature model case, shown in Fig. lb). Note, 
that this procedure also requires an assumption about 
the type of noise. 

The integral in (|^) is divergent for 7 > 1 at x = y. 
Introducing a microscopic lattice cut-off, the DST can be 
calculated |l(| to give a Dirac 5-function centered around 
t = t/2: 



for any 7 > 1 



(10) 



It is a well known result [Q that for Eq. (g) , d = d u = z 
is an upper critical dimension and separates interfaces 
which are asymptotically rough from those which are 
asymptotically smooth. Thus, the fact that for dimen- 
sions above z there is no roughening, is reflected by 
a (5-function DST, i.e., all points of the interface will 
spend exactly half of their time above the mean height. 
The persistence exponent in this case is not really de- 
fined, since the persistence probability is zero. In a 
lattice model, one would expect corrections to scaling 
to the above result, and for the persistence probability 
to decay exponentially with time. When approaching 
7« = d u /z = 1 from below, the persistence exponent 
diverges; 7 = j u being a marginal case for which no nu- 
meric or analytic results have been produced yet (on per- 
sistence properties) . It is precisely the case of the two di- 
mensional EW equation with white noise. Note that the 
EW equation in any integer dimension (d > 1) with vol- 
ume conserving noise is in the smooth phase, and would 
therefore have a 5-function DST (7=1 + d/2), with the 
persistence exponent undefined (or formally infinite). 
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FIG. 2. The behavior of DST as a function of 7 for eq. 
([]), with white noise. Each inset shows the function F 1 {x) 
vs. x. For 7 = 0.25 and 7 = 0.5 we simulated the Langevin 
equation with white noise in one dimension at z = 4 and 
z = 2, respectively. For z = 2 the simulation parameters were 
the same as for the thick line in fig. la) and for z = 4 they 
were the same as for the continuous line in fig. lb) 

From the scaling relation (Q) one can infer the exis- 
tence of a new critical 'dimension' 7* both for the de- 
terministic and noisy cases: since the tails of DST give 
the persistence probability, which has a power law decay 
(P ± ^ t~ e ), the scaling function F 1 must obey the be- 
havior F 1 (x) ~ x 6 " 1 , for iCl, and 1 — x <C 1 in order 
that (0) be satisfied. For 9 < 1 the DST has integrably 
divergent tails while for 9 > 1 the tails vanish (as x " 1 ). 
In the former case the sites are more likely to be found in 
a positive or negative persistent state (i.e. with a height 
that did not change sign at all) , while in the latter case 
persistent sites will be an absolute minority (with van- 
ishing measure as t — > 00). Since f? 7 is a monotonically 
increasing function of 7 the equation 7 = 1 will be sat- 
isfied at a unique value of 7*. At this value F 7 » is flat 
at the tails, it neither falls to zero nor diverges. This 
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shows that the value of 7 . = 6 y d , = 1 is special. It is 
possible that F 7 * can still have some structure around 
r/t = 1/2, but the simplest possibility is that it is a top- 
hat function. In this case \ii = 1/3, and therefore 7* can 
be calculated after (numerically) inverting G(^f) = 1/6, 
using Eqs. (H) and (^). For the deterministic case one 
obtains 7* = 17.983.. and for the noisy case 7* = 0.438... 
For the noisy case, our numerical simulations are com- 
patible with 0.25 < 7* < 0.5, as can be seen from Fig. 2. 
It is interesting to note that the permanent presence of 
noise 'brings down' this critical 7* to a sub-unitary value 
as compared to the deterministic case. 

The DST for 7 = is exactly known, and is called the 
'arcsine law' in the mathematical literature jy]: Fq(x) = 
l/(n-\/x(l — x)), which can also be derived from Eqs. 

-(0) presenting a novel alternative to this venerable old 
problem. Fig. 2 summarizes our findings on the different 
regimes for the DST of the noisy case of Eq. (||) with the 
two critical 'dimensionalities' 7* and j u - 
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formation and coarsening. The DST mirrors all these 
morphological and structural characteristics. Nonlinear- 
ity is obvious from the asymmetric shape. The right end 
of the curve has the highest value, meaning that the sites 
are most likely to be found in a positive persistent state, 
i.e., they belong to structures that stayed above the mean 
height all the time, namely stable mounds. On the con- 
trary, the left end when compared to the right one, is in 
minority, showing that the stable crevices, or valleys will 
contain only a small fraction of the sites, which points 
to a mounded morphology with high skewness. The fact 
that a site has a small probability to survive for a long 
time in a crevice, means that the valleys tend to dis- 
appear during time-evolution, i.e., there must be coars- 
ening. This shows the intimate connection between the 
coarsening and persistence properties of a interface mor- 
phology, which is the topic of a separate, forthcoming 
publication. 

In summary, the DST proves to be very sensitive to 
the details of the morphological dynamics, and can pro- 
vide crucial information on the non-equilibrium interface 
fluctuations. 
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FIG. 3. a) DST for the d = 1 KPZ equation with white 
noise at t = 512, t — 1024 and for bimodal noise at t — 2048, 
b) DST for the SOS DT model with Schwoebel barriers at 
t = 10 3 steps shown for three different values of the parameter 
p u which is the probability for an atom to attach to a lower 
step. The system size was L = 100, and the averaging was 
made over 2 ■ 10 4 runs for each curve. 

We shall now present and briefly discuss the numeri- 
cally obtained DST's for two non- linear systems: 1) the 
one dimensional KPZ equation and 2) the Das Sarma- 
Tamborenea (DT) SOS model with Schwoebel barriers. 
Fig. 3a shows the KPZ case at different times and with 
two different noise types (Gaussian and bimodal), using 
the discretization scheme introduced in Ref . |lj] . For the 
case of Gaussian noise one can see that the DST satisfies 
the general scaling form (Q) but with an asymmetric scal- 
ing function F(x) since the h — > —h symmetry is broken, 
reflecting the nonlinear character of the KPZ equation. 
The DST for the case of bimodal noise has a different 
shape to that for Gaussian noise, which was still evolv- 
ing for the largest times we observed (t ~ 10 4 ), indicating 
either very long crossover times, or else a more compli- 
cated scaling form. Fig. 3b shows the DST obtained nu- 
merically from the DT model with Schwoebel barriers 
p3[. This system is highly non-linear, exhibits mound- 
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